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We report a neutron scattering study on the tetragonal compound Sr2Cu304Cl2, which has 
two-dimensional (2D) interpenetrating Cuj and Cu/j subsystems, each forming a S — 1/2 
square lattice quantum Heisenberg antiferromagnet (SLQHA). The mean-field ground state is 
degenerate, since the inter-subsystem interactions are geometrically frustrated. Magnetic neu- 
tron scattering experiments show that quantum fluctuations lift the degeneracy and cause a 
2D Ising ordering of the Cuii subsystem. Due to quantum fluctuations a dramatic increase 
of the Cui out-of-plane spin-wave gap is also observed. The temperature dependence and the 
dispersion of the spin-wave energy are quantitatively explained by spin-wave calculations which 
include quantum fluctuations explicitly. The values for the nearest-neighbor superexchange 
interactions between the Cu/ and Cu// ions and between the Cu// ions are determined ex- 
perimentally to be Ji-n = — 10(2)meV and Jn = 10.5(5)meV, respectively. Due to its small 
exchange interaction, Jn, the 2D dispersion of the Cu/j SLQHA can be measured over the whole 
Brillouin zone with thermal neutrons, and a novel dispersion at the zone boundary, predicted 
by theory, is confirmed. The instantaneous magnetic correlation length of the Cun SLQHA is 
obtained up to a very high temperature, T/ Jn ~ 0.75. This result is compared with several 
theoretical predictions as well as recent experiments on the S — 1/2 SLQHA. 



I. INTRODUCTION 

Quantum magnetism has been studied for many 
decades since the advent of quantum mechanics. Most 
of the early theoretical work is based on semi-classical 
methods such as spin-wave theory. Quite remarkably, 
spin-wave theory has been successful in describing many 
physical properties of a variety of magnetic systems. De=. 
spite the fact that it is essentially a l/(zS) expansion,El 
where z is the coordination number, and thus one would 



expect it to be less accurate for a small spin quantum 
number S, spin-wave theory has been a very powerful 
tool in investigating quantum magnetism, and the "semi- 
classical" description of quantum magnets has been suf- 
ficient to understand most experimental results. In a 
seminal paper, Haldane pointed out the special signifi- 
cance of the spin quantum number in one-dimenskinal 
(ID) quantum Heisenberg antiferromagnet (QHA).H In 
his now famous conjecture, he mapped the ID QHA onto 
the quantum nonlinear a model (QNLcrM), and noted the 



fundumental difference in the ground states for half-odd- 
integer S and integer S. Specifically, the fD QHA with 
half-odd-integer S has a quasi-long-range ordered ground 
state with a gapless excitation spectrum, while that with 
integer S has a quantum disordered ground state with 
a large energy gap in the excitation spectrum. Subse- 
quent developments of quantum field theory, numerical 
simulations and experiments have confirmed Haldane's 
conjecture. 

In contrast, quantum effects in the two-dimensional 
(2D) QHA are typically less dramatic. In fact, the qual- 
itative behavior-of the 2D QHA is similar to that of 
the classical one.0 Quantum fluctuations usually manifest 
themselves through uniform renormalization of physical 
quantities, such as the staggered magnetization or the 
spin-wave velocity. However, in certain magnetic sys- 
tems, where the primary exchange couplings are highly 
frustrated, the effects of quantum fluctuation can be re- 
vealed qualitatively in the low-energy spin dynamics. As 
an example, isostructural compounds Sr2Cu304Ci2 and 
Ba2Cu304Ci2, the so-called 2342 materials, offer a dra- 
matic and clear demonstration of such quantum effects 
as "order from disorder" .□ In this paper, we describe our 
detailed neutron scattering study of the frustrated 2D 
S = 1/2 antiferromagnet Sr 2 Cu304Ci2, including exper- 
imental evidence for quantum fluctuation induced order. 
Some of the results reported here were briefly presented 
in a recent letter. El 

The discovery of high temperature superconductivity 
in 1986 has triggered much work on the magnetism in 
lamellar copper oxides. These materials contain Cu02 
planes whose 2D spin fluctuations can be modeled by the 
S = 1/2 square lattice (SL) QHA. Through a combina- 
tion of experimental, numerical, and theoretical efforts, 
a quantitati\je„ujaderstanding of the S = 1/2 SLQHA 
has emergedBtffitJ Notably, neutron scattering measure- 
ments of the instantaneous spin-spin correlation length 
of the model compound Si'2Cu02Cl2 are found to be in 
quantitative agreement with quantum Monte Carlo re- 
sults and both in turn_are well-described by analytic 
theory for the QNLerM.Q Angle resolved photoemission 
spectroscopy (ARPES) on this insulating system has also 
provided important information on the behavior of a sin- 
gle hole in a paramagnetic background, hence testing the 
applicability of the t — J model.EE3 

The structure of Sr 2 Cu304Cl2, shown in Figjl](a), is 
similar to that of S^CuG^C^. As shown in Fig. [ijb), the 
Cu02 layers are replaced by CU3O4 layers, which contain 
an additional Cujj ion at the center of evep4.second pla- 
quette of the original CU/O2 square lattice.EIrfL3 The con- 
figuration in the neighboring plane is obtained by trans- 
lating the whole plane by (| |). The in- plane isotropic 
interaction <//_// between Cuj and Cu// subsystems is 
frustrated such that they form interpenetrating S = 1/2 
SLQHA's with respective exchange interactions Ji and 
J ii- 

Due to the complete frustration of the isotropic cou- 
pling between Cu/ and Cu//, 2342 exhibits many fasci- 



• Cu 2 r + 




FIG. 1: (a) Crystal and magnetic structure of 
Sr2Cu304Cl2. Ordered spin directions for copper spins 
are shown as arrows, (b) CU3O4 plane and various ex- 
change interactions between spins. 



nating magnetic phenomena. In their magnetic suscep- 
tibility and electron paramagnetic resonance measuije=. 
ments on Ba 2 Cu304Cl2 powder, Noro and coworkersliij 
first observed anomalous features at T ~ 320K and 
T ~ 40K and attributed these to respective antiferro- 
magnetic ordering of the Cu/ and Cu// spins. Subse- 
quent neutron scattering measurements by Yamada et 
a/.ta showed that 2342 exhibits antiferromagnetic order 
of the Cuj and Cu// subsystems below the respective 
Neel temperatures Tm,i q«|d Tnji- Far-infrared elec- 
tron spin resonance £ESR)t3 and submillimeter wave res- 
onance experimentsc3 showed that there is a low energy 
out-of-plane excitation in the long- wavelength limit. The 
dispersion of a single hole in both antiferromagnetic and 
paramagnetic spin background was measured simulta- 
neously in the . same .CU3O4 plane by ARPES experi- 
ments on 2342.IIZIBII3 One of the most intriguing fea- 
tures in earlier studies is the |yegki ferromagnetic mo- 
ment that appears below T/v^MoEjEJ We have recently 
reported that anisotropic bond-dependent interactions 
such as pseudo-dipolar coupliiigjj-^an in fact explain such 
weak ferromagnetic momcnts.l 



We consider two specific consequences of the frustra- 
tion in this paper. First, in the mean field approximation, 
the Cuj and Cu// subsystems are decoupled, so that in 
addition to the well known Cu/ SLQHA, the Cu//'s form 
their own S = 1/2 SLQHA with an order_pf magnitude 
smaller superexchange, Ju. Chou et alE!i have shown 
that the magnetic susceptibility of the Cu// subsystem 
is very well described as a S = 1/2 SLQHA by compar- 
ing the experimental result with the results of a quantum 
Monte Carlo calculation. Since Ju ~ 10 meV is matched 
well with the energy of thermal neutrons, this is an ideal 
S = 1/2 SLQHA system for neutron scattering experi- 
ments. In Sec. IV D, we show the spin- wave dispersion 



of the Cu/7 subsystem throughout the entire Brillouin 
zone, including a theoreticall y pre dicted dispersion along 
the zone boundary. In Sec. VB, the correlation length 
data measured from the Cu// SLQHA are presented as 
a function of temperature and compared with various 
theoretical predictions as well as quantum Monte Carlo 
results. Because Ju is an order of magnitude smaller 
than J/, we are able to access a rather high temperature 
{T/J H » 0.75). 

Second, because of the frustration, we can observe the 
direct effect of quantum fluctuations. When a system can 
be separated into two Heisenberg antiferromagnetic sub- 
lattices, so that the molecular field of the spins in each 
sublattice vanishes on the spins of the other, then within 
mean field theory the ground state has a degeneracy 
with respect to the relative orientation of the sublattices, 
and the excitation spectrum contains two distinct sets of 
zero energy (Goldstone) modes, reflecting the fact that 
these subsystems can be rotated independently without 
cost in energy. n This degeneracy is removed by fluctua- 
tions. ShenderQ showed that quantum spin- wave inter- 
actions prefer collinearity of the spins in the two sublat- 
tices. This has the following experimental consequences 
in Sr2Cu304Cl2: The symmetry of the critical fluctua- 
tions of the Cu// system is lowered to Ising due to the 
fluctuation-driven collinearity, and the spin-wave mode 
corresponding to the relative rotation of sublattice mo- 
ments develops a gap. Indeed, such a gap was considered 
in the garnet Fe2Ca3(Ge04)3. However, since a similar 
gap could also arise from crystalline magnet ie j. anis ptropy, 
the final identification was rather complex.c3nJc3 

Our inelastic neutron data in Sec. |\J show a dramatic 
increase of the Cuj "out-of-plane" gap below Tjyjj (see 
Fig. g) , which clearly reflects a coupling between the Cuj 
and Cu// spins. However, within mean field theory this 
coupling due to frustrated interactions must vanish by 
symmetry. Accordingly, we conclude that the enhanced 
gap for T < Tff ii is due to quantum fluctuations. Heuris- 
tically, the lowering of the symmetry on the Cu/j site 
due to the ordering of the Cu/'s is sensed through the 
quantum fluctuations. This identification is corroborated 
further by detailed theoretical calculations, which use 
parameters determined independently, alheit less accu- 
rately, by the susceptibility measurements. c3 

There have been numerous studies on this peculiar or- 



def-Jrom disorder effect on various systems. Villain et 
alxn studied a generalized frustrated Ising model in two 
dimensions and found that the system does not have long 
range order at T = 0, but is ferromagnetically ordered 
at low but non-zero temperature; thermal fluctuations 
are necessary to stabilize the ordered state; thus the*; 
termed this phenomenon order from disorder. ShenderQ 
showed that quantum fluctuations can also cause order 
from disorder phenomena in frustrated magnetic sys- 
tems. In addition to thermal or quantum disorder, sub- 
stitutional disorder also causes ordering in such frus- 
trated magnetic systems. HenleycH studied order from 
substitutional disorder in a planar antiferromagnet on 
a square-lattice with a strong second nearest neighbor 
exchange and discovered that anti-collincar order is sta- 
blized by substitutional disorder, in contrast with the 
collinear ground state due to thermal or quantum disor- 
der. Chandra et alr3 investigated the Heisenberg model 
on such a lattice using analogies between quantum an- 
tiferromagnetism and superfluidity. The Heisenberg an- 
tiferromagnet on the layered body centered tetragonal 
structure, where the inter-planar coupling is fully frus- 
trated, has been also studied extensively, mainly due 
to its similarity . t o tha s tructure of high temperature 
superconductors.EaE^tilO 

The result of our neutron scattering measurements on 
single crystals of Sr2Cu304Cl2 are presented in the fol- 
lowing sections. In Sec. |n| a brief des cripti on of mag- 
netic neutron scattering is given. In Sec. Ill A , our elastic 
neutron scattering results along with the spin structure 
deduced from these results are presented. The order pa- 
rameter measurements for both the Cu/ ordering and the 
Cu/7 ordering are shown in Sec. Ill B . The spin- wave 



calculations as well as our inelastic neutron scattering 
results are discussed in Sec. IV. Section IV contains a 



large amount of data; therefore, it is divided into four 
subsections. We b riefly discuss the theoretical spin- wave 
calculation in Sec. IV A. Spin- waves at T > Tjv n, where 



the Cu j-Cu/ / interaction can be ignored, are presented 
in Sec. pVBl. The Cuj-Cujj interaction at T < T N n 



and the effect of quantum fluctuations on the spin dy- 
namics are discussed in Sec. IV C. In Sec. IV D, the 2D 
dispersion of the S = 1/2 SLQHA (Cu/j subsystem) is 
presented. We have also studied the critical behavior of 
the S = 1/2 SLQHA by measuring the static correlation 
length, which is presented in Sec. M. Finally, some un- 
resolved issues and future experiments are discussed in 
Sec. 0. 



II. MAGNETIC NEUTRON SCATTERING 

A General cross section 

In a neutron scattering experiment, the key variables 
arc the neutron energy change and the concomitant 
change in neutron wave vector. We denote the momen- 



turn and energy transfer by Q and ui, which are given 
by Q = k,j — ky and u = E^ — Ef, respectively. We 
use units in which fi = ks — 1 and the scattering vector 
Q = (^-H, Zz-K, ^-L). Throughout this paper, we use 
q to denote physically relevant momentum transfer; that 
is, the momentum transfer with respect to the reciprocal 
lattice vector G: q = Q G. 

The partial differential cross section for . spm only scat- 
tering of unpolarized neutrons is given by 33 ' 



(a) domain A 



(b) domain B 



dQ.dE 



f 



hi 



/ 2 (Q)£(<W 



Q a Qp)S a P(Q,u>), (1) 



a/3 



where Q = Q/Q, and /(Q) is the magnetic form factor, 
which is the Fourier transform of the spin-density distri- 
bution around the magnetic ion, and hence depends on 

An important feature of magnetic scattering is the di- 
rectional dependence through the geometric factor {5 a p — 

QaQp), which picks out the components of the magneti- 
zation perpendicular to the momentum transfer Q. The 
quantity ^"^(Q, w), known as the dynamic structure fac- 
tor, is the Fourier transform in both space and time of 
the spin-spin correlation function. The latter is the ther- 
mal average over the correlations between the component 
along the a-axis of a spin at the origin at time zero and 
the component along the /3-axis of a spin at site r at time 
t: 

i r°° 

S Q/3 (Q, uj) = —Y dt e '(Q- r --*) (S a (0, 0)SP(r, t)}. 

(2) 

The static structure factor is obtained from the Fourier 
transformation of the equal-time correlation function, 
and measures the instantaneous correlations between the 
spins: 
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(3) 



In principle, one can obtain the static structure factor 
by directly measuring the entire dynamical spectrum 
S(Q,u>) and doing the energy integration at each Q. 
However, in most systems this is impossible within a rea- 
sonable time scale. Fortunately, in 2D magnetic systems 
such as the lamellar copper oxides, the energy integration 
is effectively done by detecting neutrons without energy 
discrimination in a special scattering geometry. One then 
can determine the instantaneous correlation function in 



one scan.i 



B Elastic scattering cross section 

For elastic neutron scattering from collinearly ordered 
magnetic moments, the scattered intensity can be ob- 
tained from Eq. (flh: 
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FIG. 2: Two magnetic domains present in S^CuaC^C^ 
at zero magnetic field. 

where e the direction of the staggered magnetization. 
There are three factors contributing to the intensity of 
magnetic Bragg peaks: the geometric factor (1 — (Q-e) 2 ), 
the magnetic structure factor Fm(Q), and the magnetic 
form factor /(Q). 

We also need to consider magnetic domains due to 
the tetragonal symmetry of the crystal. Consider the 
two types of structure shown in Fig. |[ where only Cu 
spins are shown. In a realistic single crystal in zero mag- 
netic field, these two types of magnetic domain can be 
equally populated. These two domains give rise to dif- 
ferent magnetic reciprocal lattice vectors. As shown in 
Fig. pi different domains give different Guu magnetic 
Bragg reflections, which will prove useful in elucidating 
the spin structure of the Cu// magnetic lattice. The Cuj 
magnetic peaks, on the other hand, only occur on top 
of allowed nuclear Bragg reflections, and do not occur in 
the (H K 0) zone. 



C Inelastic scattering cross section 

In conventional spin-wave theory for a two-sublattice 
antiferromagnet, one obtains two eigenmodes. If the 
spins are ordered in the z-direction, the two modes have 
eigenvectors in the direction of x and y. For a Heisenberg 
model, these two modes are gapless Goldstone modes due 
to the continuous symmetry. However, in the presence of 
an uniaxial anisotropy, this continous symmetry is bro- 
ken, and both modes obtain energy gaps; this energy gap 
corresponds to the energy cost in rotating the spins away 
from the z-direction. For an XY anisotropy, only one 
mode has an energy gap, corresponding to the energy cost 
in rotating spins out of the xy-plane. The other mode is 
a zero-energy mode at q=0, since the continuous symme- 
try is preserved in the xy-plane. Since the polarization 
of the eigenvector of the gapped mode is perpendicular 
to the xy-plane, we call this mode an out-of-plane mode, 
while the gapless mode is called an in-plane mode. 

The direction of the eigenvectors plays an important 
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FIG. 3: Reciprocal lattice diagrams of the three different 
2D zones employed in our experiment: (a) (H K 0) zone, 
(b) (H L) zone, or equivalently, (0 K L) zone, and (c) 
(H H L) zone. Note that Cuj magnetic Bragg reflections 
coincide with nuclear Bragg positions. 



role in the neutron scattering cross section. By consid- 
ering geometric factors for both domains in Fig. g, one 
can show that the inelastic cross section from spin-waves 
reduces to 
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where </> is the angle subtended by Q and [0 1], and 
5" and S 1 - denote the dynamic structure factor of the 
in-plane and out-of-plane spin-wave modes, respectively. 
The out-of-plane component of the dynamic structure 
factor, S (Q,w), is well approximated by 
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where n(u)±) = l/(e UJ± ^ T — 1) is the Bose population 
factor, r _1 is a small magnon lifetime, and lu± is the 
out-of-plane gap. A similar relation holds for the in-plane 
component S"(Q,co) with uiu replacing ui±. 



D Experimental details 

We have carried out both inelastic and elastic neutron 
scattering experiments with the triple-axis spectrometers 



at the High Flux Beam Reactor (HFBR), Brookhaven 
National Laboratory, and at the National Institute of 
Standards and Technology, Center for Neutron Research 
(NCNR). Our measurements were done mostly on ther- 
mal beamlines at these facilities, except for the data 
shown in Fig. |l2j, which were obtained using cold neu- 
trons. Large (dimension 2 x 2 x 0.5 cm 3 ) single crystals 
of Sr^CuaO^CLj, grown by slow cooling of a melt con- 
taining a CuO flux, are used in the experiment. The 
crystals remain tetragonal (space group H/mmrri) for 
\hK < T < 550K with lattice constants a = 5.457A 
and c = 12.521 at T < 50K.E3 

The (002) reflection of pyrolytic graphite (PG) was 
used as both monochromator and analyzer. A PG filter 
was placed either before or after the sample to elimi- 
nate higher order contamination. Various experimental 
configurations with different sets of collimations and neu- 
tron energy were used. A typical setup used in the in- 
elastic experiments was a fixed final neutron energy of 
E f = 14.7 mcV and collimations of 40'-40'-Sample-40'- 
80'. The sample was sealed in an aluminum can filled 
with helium exchange gas, and mounted in a closed-cycle 
helium regrigerator. The temperature was controlled 
within ±0.2K in the range 10K < T < 400K. 



III. ANTIFERROMAGNETIC ORDERING OF 
COPPER SPINS 

Since the 2D SLQHA does not have long-range order 
at T > 0, such order must arise from spin anisotropy 
terms or inter-plane coupling. For S^CuO^C^, as the 
Neel temperature is approached from above, succes- 
sive crossovers from 2D Heisenberg to 2D XY to three- 
dimensional (3D) XY behavior are expected to take 
place, cJJpeit with a 3D critical regime that is extremely 
narrow EO For Sr2Cu304Cl2, on the other hand, the 
inter-plane coupling between Cu/ spins, Ji^d, is larger 
than the XY anisotropy. Upon lowering the temperature 
in the paramagnetic phase, we then expect a crossover 
from 2D Heisenberg behavior, characterized by a spin- 
spin correlation length ^o(^) that increases exponentially 
in T _1 , to 3D Heisenberg behavior at a temperature 
given by the relation £,qJi,3d/Ji ~ L whete Ji.zd is 
the interplane coupling between Cu/ spins .Ej In other 
words, we expect 3D effects to become important for 
£,o/a ~ 30. The correlation length of the SLQHA-irk 
known to be about 40 lattice constants at T/J ~ 0.26BE3 
For Jj ~ 130 meV, this corresponds to a temperature of 
~ 390K, which agrees with Tnj. 

Unlike for the Cu/ subsystem, the isotropic inter- 
plane Cu//-Cu« coupling is frustrated, similar to that 
of Sr2CuO2Cl2 therefore, Tnji is expected to be de- 
termined mostly by spin anisotropies, originating from 
both in-plane quantum fluctuations and from inter-plane 
dipolar and pseudodipolar interactions. For T < Tjv,j, 
the ordered Cuj spins fluctuate mainly in the directions 



transverse to their staggered moment M S] j. Ji-n then 
generates fluctuations in the Cu// spins along the same 
direction, causing an effective reduction in the.-corre- 
sponding transverse exchange components of Juu This 
yields an effective biquadratic term — S(Si ■ S//) 2 , where 
8 oc J]_ H /{Ji + Jn). This implies an effective Ising-like 

anisotropy Juttfj oc (5S 2 , which favors ordering of the 
Cu/7 spins collinearly with S / , consistent with our mea- 
sured structure, Fig. 0(a). Indeed, Tnji ~ 40K agrees 
with £o(Tnji) 2 o;jj ~ 1, where ajf ~ 0.01 is indepen- 
dent ly ded uced from our spin-wave gaps, as discussed in 
Sec. |TVD . We next show experimentally that the order- 
ing direction of Cu// spins is indeed parallel to that of 
the Cu/ subsystem, and that this ordering is a 2D Ising 
transition. 
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A Magnetic structure 

The Cu/ spin ordering direction shown in Fig. 
|l| has been .determined in previous magnetization 
measurements£3 This ordering is similar to that in a 
bilayer cuprate Yl^CusOg- Unlike other "214" type 
materials, the Cu/ spins of "2342" have unique nearest 
neighbors in the c-direction, just like YBa 2 Cu 3 06. For 
such a structure, the observed ordering direction along 
the Cu-O-Cu bonds has been attributed to the quan- 
tum fluctuations .c2l 

Next, let us consider Cu// magnetic diffraction peaks 
in the (H K 0) zone. In Table 1, we summarize our re- 
sults. The neutron energy was fixed at 14.7 meV, and 
collimations of 20'-40'-S-40'-80' were used. One can fit 
the peak intensities with Eq. (0) with e as a free param- 
eter. The fit gives the spin direction shown in Fig. @; 
namely, e along the [1 1 0] direction for domain A, and 
along the [110] direction for domain B. In order to de- 
termine how the copper oxide layers are stacked, we show 
in Fig. EJthe peak intensity for each magnetic Bragg peak 
in the (H H L) zone. From the domain structure and the 
stacking scheme in Fig. |2|, one can show that magnetic 
Bragg peaks occur at even L due to domain A, together 
with magnetic Bragg peaks at odd L due to domain B. 
Since our momentum transfer, Q, is along the [H H L] 
direction, Q is always perpendicular to the spin ordering 
direction, e, in domain A. However, in domain B, this 
is not true, and e • Q = sin 0. Therefore, the geomet- 
ric factor only matters for the peaks from domain B. As 
shown in Fig. EL only the odd-L data exhibit the expected 
geometric factor dependence. In fact, the agreement is 
excellent between odd-L data (triangles) and the calcu- 
lation (x). It should be noted that the stacking scheme 
of the Cu// spins is identical to that of Si^CuC^Ck-Jid 

In their study on YBa2Cu306.i5, Shamoto et alc3 re- 
ported an anisotropic Cu magnetic form factor, which 
depends not only on the magnitude of Q, but also on 
the direction of Q. Specifically, the magnetic form factor 
was found to drop more rapidly with increasing Q, if Q 



FIG. 4: Cu// magnetic Bragg peak intensity for (H H L) 
is plotted for even-L (open circle) and odd-L (closed tri- 
angle). The solid riipe is a plot of magnetic form factor 
of free Cu 2+ ions.L3 The symbol x represents the mag- 
netic form factor squared multiplied by the geometric 
factor in Eq. (Q). Clearly, even-L data show the same 
Q-dependence as the magnetic form factor, while odd-L 
data do not. The agreement between odd-L data and the 
calculated results are very good. 



is perpendicular to the L-direction. The small deviation 
between the even-L data (open circles) and the solid line 
at large Q is probably due to this anisotropy in the mag- 
netic form factor, since most of our large-Q data has a 
relatively small L-component. 

From the magnetic Bragg peak intensity, we have 
also estimated the value of staggered magnetization at 
T = 10K as M sJ « 0.4(2)/x B and M sJI « 0.8(2)^ B for 
Cu/ and Cu//, respectively. Within experimental error 
bars these values reasonably agree with the theoretically 
expected value 0.6^_b.e3 



B Order parameters 

The antiferromagnetic Bragg intensity is proportional 
to the square of the staggered magnetization, M s , which 
is the order parameter of the Neel ordered phase. We 
measured the temperature dependence of the (1 1) 
peak, using neutrons with energy 13.4meV and collima- 
tions of 20'-40'-S-40'-80'. The temperature dependence 
of the magnetic Bragg intensity at the (10 1) reciprocal 
lattice position is shown in Fig. ||(a). Since nuclear Bragg 
scattering is only weakly temperature dependent, we sub- 
tract the high-temperature (10 1) nuclear intensity from 
the observed intensity. We also studied the temperature 
dependence of the (3 1) peak, which shows the same 
temperature dependence as the (10 1) peak. 

We fit the T > 300K data to the form I~(T N - T) 2/3 , 
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FIG. 5: (a) Filled circles are the integrated intensity of 
the 3D magnetic Bragg peak at (1 1). Open circles are 
the intensity on the 2D rod of Cu/ at (1 0.55). (b) 
The peak intensity of (1/2 1/2 0) peak is plotted as filled 
circles. Open circles are the intensity on the 2D rod of 
Cu/7 at (1/2 1/2 0.25). Solid lines are fits to - {T N -T) 213 
with the fitting parameters shown. 

TABTE I: Cu// magnetic Bragg intensity in (H K 0) 
zone at T — 12K. The observed intensity, I(obs), is neu- 
tron counts per minute. The intensity, I(cal), is cal- 
culated assuming the spin structure shown in Fig. @. 
The crystallographic reliability index is R = J2 |VX>bs — 

V / 7cal|/EVTobs = 0.12. 
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and obtained T NJ = 386(2)K and Br = 0.28(3); this re- 
sult is plotted as a solkLline in Fig. ||(a). 3i is consistent 
with that of LaaCu04,E3 while somewhat larger than that 
of Sr 2 Cu0 2 Cl 2 llWe associate /?/ with the 3D XY univer- 
sality class, since the inter-plane coupling drives the 3D 
ordering in Sr2Cu304Cl2- On the other hand, the domi- 
nant spin anisotropy driving the ordering in Si"2Cu02Cl2 
is the XY anisotropy within the plane; thus the system 
is presumably closer ta-tJre 2D XY universality class at 
least not too near TV .0a Ref. 45 also suggests that such 
XY systems might exhibit an order parameter exponent 
B close to \ . Another possible explanation follows experi- 
ments on La2Cu04, which were interpreted as being near 
a tricritical point, which would have exactly 3 — ^e3. 
However, it is not yet clear why many quasi-2D Heisen- 
berg and XY systems happen to be near tricriticality. 

For the Cu/j order parameter, we measure the inten- 
sity of the (| | 0) reflection as a function of temper- 
ature. The neutron energy is fixed at 14.7 meV, and 
10'-40'-S-40'-80' collimations are used. The peak inten- 
sities are shown in Fig. g(b). As in the Cu/ case, the 
solid line is the result of fitting data for T > 30K with 
T N . H = 39.6(4)K and Bxx = 0.13(1). This strongly sug- 
gests that the ordering of the Cu// spins is in the 2D Ising 
universality class. As discussed in Sec. |, this provides 
strong evidence of fluctuation driven ordering. However, 
one cannot, from this measurement alone, rule out the 
possibility of the presence of crystalline Ising anisotropics 
in the Cujj-Cujj superexchange interactions. As shown 
in Sec. IV, evidence from spin dynamics experiments is 



necessary to clarify this point. 

Another way to probe the 3D magnetic ordering is to 
plot the intensity of the 2D magnetic rod as a function 
of temperature. Quasi-2D materials, such as K2N1F4, 
La2Cu04,c3 and Sr 2 Cu02Cl2j3 show strong 2D dynamic 
fluctuations above the 3D ordering temperature Tjv; this 
is exhibited as rods of scattering perpendicular to the 
2D plane, whose locations are shown in Fig. 0. One can 
observe this by accepting all energies of neutrons at the 
detector in the two-axis configuration. At TV, the 2D 
inelastic scattering intensity begins to decrease rapidly, 
as the spectral weight is shifted from 2D inelastic scatter- 
ing to 3D Bragg scattering. Therefore, this measurement 
shows the 2D nature of the system, as well as the 3D tran- 
sition temperature, complementing the order parameter 
measurement. In Fig. @(a) and (b), we show these 2D rod 
intensities at (1 0.55) and at (0.5 0.5 0.25) for the Cu/ 
and Cu// sublattices, respectively. Indeed, we see rapid 
decreases of both the Cu/ and the Cu// 2D rod intensi- 
ties as the system is cooled through their respective 3D 
Neel transitions. The non-zero intensity below the Neel 
temperature is due to the contributions from phonons. 
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TABLE II: Parameters used in the spin Hamiltonian [Eq. 
(0)]. These values are determined from our neutron scat- 
tering experiment. Superexchange energies are in units 
of meV and a is dimensionless. 



meaning 

Ji Cuj-Cu/ superexchange (in-plane) 

Ju Cuj_f-Cu/_f superexchange 

Ji-ii Cui-Cun superexchange 

Jl,3D Cui-Cui superexchange (out-of-plane) 

ati XY-anisotropy in Jj (T=200K) 

an XY-anisotropy in J u (T=10K) 



value 
130(5) 
10.5(5) 
-10(2) 

0.14(2) 

5.2(9) x 10~ 4 

1(5) x 10~ 4 



IV. SPIN DYNAMICS 

The spin Hamiltonian used in the spin-wave calcula- 
tion and the data analysis is as follows: 

H = Hi + Tin + 7~t mt (7) 

Hi = Ji 2_^ (Si • Sj - aiSfSj) + Ji^d 2_^ S, • Sj 



: ;j)l 



(i,j)l,30 



Till — Jll 2^, (^m ' ^ n ~ a H^m^n) 
(m,n)ii 

Tiint = Jl-II 2_^ Si ' ^ m ' 

where i,j and m,n denote Cu/ sites and Cu// sites, re- 
spectively. The symbols (i,j)j and (i,j)i^D label Cu/ 
intra-planar and inter-planar nearest neighbors, whereas 
(m,n)// and (i,m)i-n refer to the nearest-neighbor 
Cuij-Cujj and Cu/- Cu/j bonds, respectively. The re- 
duced exchange anisotropy, a — (J— J z )/J, is used here, 
and is therefore dimensionless. We left out other smaller 
terms, such as the in-plane anisotropy in Jj and Ju, 
the pseudo-dipolar interaction between Cu/ and Cu/j, 
J p d, the interplanar dipolar Cu/j .-Cujj interaction, and 
the four-fold anisotropy term.EJEa It turns out that these 
small terms do not affect the spin dynamics on the en- 
ergy scale probed by thermal neutrons, although they are 



obtained from our data analysis, the details of which are 
discussed next, are summarized in Table |l|. 



A Spin-wave theory 

1 T=0 

Our measured spin-wave energies can be explained 
within the framework of T — interacting spin-wave 
theory (SWT). The theory is discussed in detail in Ref. 
|49| . Here we give only a brief summary with the salient 
results. Starting from the spin structure shown in Fig. 
|l|, we express each of the six spins (four Cu/'s and two 
Cu// 's) in the unit cell by the Dyson-Maleev transforma- 
tion for general spin S. The sums in Tii and Tin are then 
truncated at the harmonic order in the spin-wave boson 
operators. However, the Tii n t term vanishes at the zone 
center, and therefore has effects only if one expands it to 
quartic order. We then approximate each product of four 
spin-wave operators by contracting operator-pairs in all 
possible ways. This yields new quadratic terms, whose 
coefficients contain the parameter d = 2Ji-u(ae)/S, 
where a and e are boson operators associated with Cu/ 
and Cu//, respectively. This coefficient contains the fac- 
tor 1/S, thus representing quantum corrections due to 
spin-wave interactions. The spin-wave energies are then 
found as the eigenvalues of the 6x6 matrix which -arises 
from the resulting bilinear spin-wave Hamiltonian.E3 

Since the magnetic unit cell contains 6 Cu spins, the 
spin-wave spectrum has six branches. Two of these are 
optical modes which are practically degenerate at lu — 
ASZ c Ji. In this paper we will only discuss the remaining 
four modes. The large in-plane spin-wave velocity for 
the Cu/ spins makes it difficult to study the dispersion 
other than at the 2D zone center along the L direction, 
where the mode energies can be found analytically. The 
energies of these modes at T = for wavevectors (1 L) 
are (in order of increasing energy) 



essential in explaining such behavior, as, 
sition or ESR experiment results. E2lr3l 



re spin-flop tran- 
The parameters 



u>i 



--S^/32J II Sx 3 /(S + 2x 3 ) 



2JiZ 2 ai + x 3 



. 2 =S x lmJu[Zp u a II + 8- Aj iZ2 y 



2x-. 



u> 3 =S\ 8Ji 



uj4 =S\ 8Ji 



2x 3 + 6 [1 



Ji~ii Ju 25 
J i J i S + 2x3 



AZ 2 g J I a I + 2x ?i + 5[ 1 



Ji-ii _j_ Jjj_ 
Ji 



2d 



J/ 6 + 2x 3 + AZlJ T ai 



(8) 

(9) 

(10) 
(11) 



where Z g = 1 + 0(1/5*) sa 0.6 is the quantum renor- S — 1/2, and x 3 = ZfJ/^fl + cos(-7rL)], where Z 3 = 
malization factor for the spin-wave anisotropy gap when 1 + 0(1/S) ~ 0.9. In Eqs. (|8|-|ll|) we have kept only 



terms up to 0(1/ S). Since 6 = 0(1/ S), this term is not 
renormalized. Note that the dispersion of uii and 0J2 is 
of order 6, and hence is purely fluctuational. Note also 
that ai and Ji^d appear always with the renormaliza- 
tion factors Z g and Z3; thus, we can only determine the 
products ZgCti and Z 2 Ji^r>- 

The physics of these modes can be deduced from the 
structure of the mode energies. Because only 0J2 and 
W4 involve the XY anisotropics, ai and a//, we see that 
these modes are out-of-plane modes, i. e. modes in which 
the spins oscillate out of the easy plane. Correspondingly, 
<x>i and L03 are in-plane modes in which the spins move 
within the easy plane. In this connection note that the 
energies for u>i and u>3 can be obtained from 102 and 0*4, 
respectively, by omitting all factors which involve the XY 
anisotropics. Likewise, modes u>3 and 0J4 involve Ji and 
are hence modes which primarily exist on the Cu/ sublat- 
tice (Cu/ modes), whereas modes u>i and L02 involve Ju 
and are modes which primarily exist on the Cu// sublat- 
tice (Cu/7 modes). From this it follows that the modes 
u>i and 102 w ih have high intensity near Cu/j Bragg po- 
sitions and low intensity near Cuj Bragg positions and 
conversely for the modes W3 and 0J4. Finally, we should 
point out that u\ — > as q — > only because we have 
here neglected the small pseudodipolar interactions and 
four-fold anisotropy which lead to in-plane anisotropy. 

One should be careful in determining the absolute 
value of the XY-anisotropy of the exchange coupling, 
since the quantum renormalization factor for the spin- 
wave gap (Zg) is different from that of the spin- wave 
velocity (Z c ). Moreover, the value of Z g is not known 
accurately. Z g was first discussed by Barnes et al£j 
in their Monte Carlo study of a Heisenberg-Ising anti- 
ferromagnet. They discovered that the anisotropy gap 
was almost a factor of two smaller than that of the lin- 
ear spin-wave prediction. In their series expansion study 
of the Heisenberg-Ising model, Zheag et alE3 calculated 



Z g = 0.635(10); Singh and GelfanoS also calculated the 
renormalization of the Ising gap using the series expan- 



sion method: Z g w 0.56. This value agrees with Monte 
Carlo data from Ref. ^(J. On thq-athcr hand, Z c has 
been known since Oguchi's work,E3 and a number of 
high-precision calculations of Z c have become.-available 
recently .-The series expansion results by Singha and by 
Igarashio are Z c = 1.176 and Z c = 1.1794, respectively. 
We use the Monte Carlo result of Beaxd et at, Z c = 1.17B 
As pointed out by Barnes et a/.,t2l the factor Z g can 
be physically, understood by considering the effect of 
fluctuations.^ Spin wave theory assumes a classical Neel 
ground state with a perfectly ordered moment, and the 
resulting dispersion relation is for the spin-waves propa- 
gating in such a background. Clearly, both quantum fluc- 
tuations and thermal fluctuations substantially reduce 
the ground-state alignment; long wavelength spin- waves 
thus see a "softened" antiferromagnetic background, and 
we see the renormalization of the gap, which is propor- 
tional to the staggered magnetization that is reduced 
from its classical value due to fluctuations. The effect 
of thermal fluctuations is well known from the studies 
of K 2 NiF4 and Sr 2 Cu02Cl2, where the gap energy fol- 
lows the temperature dependence of the order parame- 
ter. Analogously, quantum fluctuations also reduce the 
gap energy from the classical value even at zero temper- 
ature. In the S = 1/2 SLQHA, the zero temperature 
ordered moment is reduced by ~ 40%. We use the renor- 
malization Z„ ~ 0.6 from this fact. 



2 Temperature dependence of the Mode Energies 

As we shall see, fitting the experimentally determined 
mode energies to the expressions of Eqs. (p[]Tl]) suggests 
that the temperature dependence of 6 is the same as that 
of M% u . Combining the zero temperature results with 
the random phase approximation results for 5 = we 
propose to describe the mode energies at nonzero tem- 
perature (but for X3 = 0) by 



Wo = Z2JhSii ZgJjiauSji 



oo^ n 



2J I Z 2 g a I S I 
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11 1 Ji-ii , Ju 



26 S 2 
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Ji J 7 5oS 2 j+AJ iai Z 2 S 2 



(12) 

(13) 
(14) 



where Sj = 5(1 - T/T NJ f', S H = S(l - T/T N . H f' 
and $o is the value of 8 for T = 0. 



ferromagnet. We also require that T <§C Tn,i, so that we 
can ignore the T-dependence of the Cuj moment. Set- 
ting 6 = 0, luj_ = L04 near the (1 L) magnetic reciprocal 



B T > T P 



At high temperatures (T 3> Tnji), we can ignore the 
Cuj/'s and treat the Cuj system as a two-sublattice anti- 
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position is given by 

lo ± = AS J! 



2Z 2 g ai- 



:(<72do) 2 



+zr^(i 



COS (7rl/)) 



1/2 



(15) 



Here, 52D is the momentum transfer in the plane, q^n = 
— y/(H — l) 2 + K 2 , and a is the lattice constant. Note 
that the distance between the Ciij— Cuj nearest neighbors 
is a/>/2. The in-plane mode, con — W3, has the same 
dispersion relation, with a.j replaced by zero, because we 
assume zero in-plane anisotropy. 

From a measurement of the spin- wave dispersion along 
[1 L], it is then possible to extract both a] and 
J 1. 3D! using Eq. rtlq). In Fig. H, such measurements 
are shown. We have chosen a relatively high tempera- 
ture of T = 200-ftf rs Tjf j/2 in order to take advantage 
of the magnon population factor n{ui) in the cross sec- 
tion. The experiment was carried out in the constant-Q 
mode, in which the final neutron energy was fixed at 
Ef = 14.7meV and the spectrometer was set to oper- 
ate in the neutron energy loss configuration. A horizon- 
tal collimation sequence 20'-40'-S-80'-80' was employed, 
which resulted in energy resolutions (full width) between 
~ 1.4 and ~ 1.7meV for energy transfers between 3 and 
12meV. In this figure, (1 1) is the zone center and (1 
2) is the zone boundary. In order to show that the gap 
of ~ 5meV at the zone center (1 1) is indeed an out- 
of-plane mode, we compared the (10 1) data with the 
scan at the (1 5), where the intensity of an out-of-plane 
mode would be reduced due to the geometric factor sin 2 cf> 
in Eq. (|. 

We have analyzed our data by convolving the spin- 
wave cross section, Eq. (|j), with the resolution function 
of the spectrometer. Attempts to fit the data with a 
single peak are not successful. Figure 0(a) shows a sum- 
mary of our results for the dispersion of the Cuj modes 
tt>3 and 0J4 at T = 200K, along the L-direction. The filled 
circles denote the in-plane mode and the open circles the 
out-of-plane mode. Fixing Jj = 130meV and S — 0, we 
fit the gapless data to W3 [Eq. (|l0|)1, and obtain Ji$d = 
0.14(2)meV. Compared to Jj ~ 130meV, the inter-planar 
coupling is rather small: Ji^d/ Ji ~ l.lxlO -3 . However, 
this Ji^zd/ Ji value is larger than those in S^CuC^C^ 
and La2Cu04, where copper oxide planes are stacked dif- 
ferently. Namely, the interplanar interaction is not frus- 
trated in Sr2Cu3 04Cl2, whereas in the other compounds 
it is frustrated. 

We then fit the out-of-plane data to Eq. ( pd| ) and find 
ai = 5.2(9) -x 10~ 4 . In their study of Sr 2 Cu0 2 Cl2, 
Greven et alu determined the XY-anisotropy in Jj as 
«2i22 = 1-4(1) x 10~ 4 , without recognizing the differ- 
ent quantum renormalization factor for the spin-wave 
gap (Z g ). If «2i22 is multiplied by (Z c /Z g ) 2 , it becomes 
5.3(4) x 10~ 4 , which is consistent with our ai. This XY 
anisotropy has been explained by Yildirim et alxl as re- 
sulting from a combination of spin-orbit and Coulomb 
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FIG. 6: Representative inelastic neutron scattering inten- 
sity from Cu/ spin waves (0*3 and LO4). Each panel shows 
a scan at T = 200K with the momentum transfer Q fixed 
as noted. The solid line is the cross section, Eq. (0), con- 
voluted with the experimental resolution function. The 
dashed lines indicate the individual contributions of the 
two spin-wave modes to the overall intensity. 

exchange interactions. 

Note that the solid lines in Fig. are produced with 
fixed Jj, a/, and J/.3_d, while the dashed lines indicate 
the individual contributions of the two spin-wave modes 
to the overall intensity. The Q-dependence of the scat- 
tering intensity only comes from the geometric factor and 
the spin- wave dispersion; thus, we are able to use a single 
set of parameters in Eqs. (Hjg) to explain the observed 
cross sections for various Q's in Fig. pi 

We next discuss the temperature dependence of the 
out-of-plane gap. In prewous studies of the tetragon. 
SLQHA K 2 NiF 4 (S = 1JB and Sr 2 Cu0 2 Cl 2 (S = l/2),l 
it has been found that the gap energies exhibit the same 
temperature dependences as the respective order param- 
eters throug hout t he entire ordered phase. As discussed 
above in Sec. [V A, this is due to the softening of the anti- 
ferromagnetic background by thermal fluctuations. Fig- 
ure shows constant-Q scans at various temperatures at 
the Brillouin zone center (1 1). The data can be fitted 
well with Eq. (0); the fitting results are shown as solid 
lines in the figure. Due to the steep in-plane dispersion, 
the fitted gap energy, indicated by the arrows, is slightly 
smaller than the apparent peak position, and the peaks 
appear to have wider widths than the resolution. 

The summary plot in Fig. M clearly shows the change 
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FIG. 7: (a) Spin-wave dispersion of lu^ and L04 along 
(1 L) at T=200K. Filled circles denote 0*3 and open 
circles L04 extracted from fitting the data. The solid lines 
are the dispersion relations, Eq. ( |l0| ) and Eq. (O). (b) 
tt>3 and 0J4 at T=30K, and 0J2 a t T—12K. The solid lines 
are the dispersion relations, Eq. (j8|-|ri|). u>2 is shown in 
open triangles. Filled diamonds are results from antifcpj 
romagnetic resonance experiment by Katsumata et alcB 

of the gap energy as a function of temperature. As ex- 
pected, UI4 follows the Cu/ order parameter from T^ji ~ 
40K up to Tnj. However, a large increase is observed be- 
low Tnjj. Further new low-energy features appear be- 
low 40K. We discuss this dramatic behavior of the long- 
wavelength spin-waves below 40K in the next subsection. 



C T < T N ,n 

The dramatic behavior of the low-energy, long- 
wavelength spin waves at T < Tjvj/ results from the 
effective biquadratic interaction, S, produced by quan- 
tum fluctuations. We now discuss the effect of S on the 
spin- wave energies on the basis of Eqs. (|[0). The en- 
ergy of the out-of-plane mode U4 increases dramatically 
as the Cu/7 spins order and S comes into play. For the in- 
plane mode the effect of S is even more dramatic because 
in the absence of in-plane anisotropy, its energy is zero for 
5 = 0. The existence of the out-of-plane mode 0J2 requires 
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FIG. 8: Constant-Q scans at the 3D zone center (1 
1) at various temperatures. Solid lines are fits to Eq. 
(||). Arrows denote the gap energies extracted from the 
fits. The final neutron energy was fixed at 14.7 meV 
and collimations of 20'-40'-S-20'-40' were used. At high 
temperatures, tighter collimations were used to improve 
the energy resolution. 



long-range order of the Cu// subsystem. Here nonzero 8 
causes an increase in the out-of-plane anisotropy (from 
Jiiotn to approximately Juan + 2Jjai) because the 
quantum fluctuations strongly couple the two subsys- 
tems. As mentioned above, the effective biquadratic ex- 
change does not create an energy gap in the mode uj\ . 



1 The Cui-like modes 

The peak in our data for T < 40K, in Fig. 0, is iden- 
tified as an overlap of peaks from the W3 and L04 modes. 
These spin- wave modes could not be resolved due to both 
the steep in-plane dispersion of the Cu/-like mode and 
the existence of a nearby phonon peak. However, one 
can obtain indirect evidence for the correctness of this 
description by exploiting the different polarizations of the 
0J3 and LU4 modes. In Fig. [LQ(a), we compare scans at the 
(1 I) position and the (10 7) position. The data at the 
(10 1) position can be satisfactorily fitted with both a 
single peak and two peaks. Specifically, the dashed line 
in Fig. |l^(a) assumes S is zero in Eq. (p^|), so that there 
is only one energy gap from L04. The solid line assumes 
non-zero 5, thus producing a double peak feature: both 
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FIG. 9: Temperature dependence of the spin-wave gap 
at the 3D zone center (1 1). Open triangles, filled 
circles and open circles denote W2, u>3, and W4, respec- 
tively. Solid, dashed, and dot-dashed lines represent re- 
spective spin- wave calculation, Eq. ( |12|JL4[ ). Inset: Same 
data plotted in a different scale to magnify the low tem- 
perature region. 



W3 and 104. Using the same set of parameters obtained 
from fitting the (10 1) data, we plot the solid and dashed 
lines for the peak profile at (1 7). Evidently, an in-plane 
gap (0J3) is necessary to explain the data at (1 7), where 
the contribution of the out-of-plane mode becomes very 
small due to the geometric factor. Therefore, we have 
shown that this peak below 40K results from an over- 
lap of u>3 and 0J4. We have, therefore, fitted all of the 
data assuming that there are two modes. We emphasize 
that for T < T^ji, the non-zero energy of ^3(q = 0) 
is a pure quantum effect; the close values of 103 and W4 
simply reflect the fact that the effective anisotropy asso- 
ciated with 5 is larger than the intrinsic Cuj out-of-plane 
anisotropy aj [see Eqs. ( |l0| ) and (|TT|)] . thus illustrating 
the quantitative importance of quantum fluctuations. 

We have measured the dispersion along the L-direction 
of the UI3 and W4 modes for T < Tjyji. Each scan is fitted 
with the cross section containing both W3 and L04. The 
fitting results are shown in Fig. ffl(b) as filled and open 
circles for 103 and u>\, respectively. The solid lines in Fig. 
0(b) are drawn using Eqs. (|^|l|) and (^||), with S = 



FIG. 10: (a) Constant-Q scan of the spin- wave gap at the 
3D zone center at 10K. The (10 1) scan shows an overlap 
of in-plane (1J3) and out-of-plane {0J4) mode, while the 
scan at (1 7) is almost entirely in-plane mode, due to 
the geometric factor of the neutron cross section. The 
baseline of the (1 7) data is offset by 50. The solid 
and dashed curves are fits for two peaks and one peak, 
respectively, (b) Dispersion along L of Cu/-like modes 
at T=30K. 



0.26(4)meV determined by fitting 0J3 with fixed J] — 
130meV. Note that we have assume d the t emperature 
dependence of S as discussed in Sec. [IV A 2 . Using the 
theoretical relation So = 0.3372J|_ // /Jj = 0.26 from 
Ref. H|, we obtain |J/_jj| = 10(2) meV, in excellent 
agreement with the earlier magnetization study.E3 



2 The Cun-like modes at (1 1) 

The low-energy mode that appears at temperatures 
below 40K is attributed to 0J2 ■ At least two experimental 
observations support this identification. In Fig. ffl(a), 
we compare this mode at different L positions; the peak 
evident at the (1 1) position disappears at the (1 7) 
position, thus proving that this gap is an out-of-plane 
mode. Next, in order to show that this mode is Cu //- 
like, the (1 1) scan is compared with the scan at the 
(1.02 1) position in Fig. |TT|(b). Although there is no 
discernible peak at (1.02 1), the remaining intensity 
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is consistent with the calculation using the Cu// spin- 
wave velocity (~95meVA). If the Cuj spin- wave velocity 
(~830meVA) is used instead, the dashed line is obtained, 
which is basically at the background level. Therefore, this 
low-energy feature is the Cu//-like out-of-plane mode: 

L02- 

In Fig. [y](c), the temperature dependence of the u>2 
gap at the (1 1) position is shown. As expected, the 0J2 
gap vanishes for T > 40K. The temperature dependence 
of ui-2, is summarized in Fig. as open triangles. The 
lines in Fig. correspond to Eqs. (Ol-O). The agreement 
between the calculation and the experimental results over 
the entire temperature range is excellent, if one takes 
into account the inherent difficulty in resolving 074 and 
consequent large error bars for 0*4. 

The dispersion of 102 along L at T = 12K is shown in 
Fig. |ll|(d); the summary is plotted in Fig. 0(b) as open 
diamonds. The solid lines in Fig. H(b) for u>\ and L02 have 
no adjustable parameters; all the parameters have been 
determined independently from separate measurements. 
We set an = 0, since our least square fit of the data to 
Eq. (||) yields an = 0.0001(5), which is indistinguishable 
from zero. 

The uj\ mode could not be identified as a distinct mode 
in our experiment, due to the presence of an acous- 
tic phonon. Note that the (1 1) position is a nu- 
clear Bragg position as well as a magnetic zone center. 
However, in a recent study, Katsumata et aZo reported 
an observation of antiferromagnetic resonance modes at 
T=1.5K using the ESR techniques. They showed that 
there are two modes: an out-of-plane mode at 422.5 
GHz(^ 1.75meV) in good agreement with our u>2 value 
\u>2{T — ► 0) ~ 1.72meV], and an in-plane mode at 36.1 
GHz (~ 0.15meV), which is too small to be observed with 
thermal neutrons. These results are plotted as filled dia- 
monds in Fig. R 

Therefore, the combined inelastic neutron scattering 
results and spin-wave calculations in Figs. R and H clearly 
demonstrate the success of our model Hamiltonian in ex- 
plaining the observed temperature and momentum de- 
pendences of the spin-waves. 
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FIG. 11: (a) Comparison of 0J2 gap at (1 1) and (1 7), 
showing that 0J2 is an out-of-plane mode. Experimental 
configuration: Ei — 14.7 meV, collimation sequence of 
40'-40'-S-40'-80', at T = 35K. (b) Comparison of uj 2 gap 
at two slightly different in-plane wave vectors. Experi- 
mental configuration: E f = 14.7 meV, 60'-60'-S-60'- 
40', at T = 10K. (c) Temperature dependence of u>2 gap 
at (10 1). Experimental configuration: Ei = 14.7 meV, 
40'-40'-S-20'-80'. (d) Scans showing a small dispersion 
of LO2 along L. Experimental configuration: Ei = 13.7 
meV, 20'-40'-S-20'-40', at T = 12K. All solid lines are 
fits to Eq. (H). Dashed curve in (b) is the same as that 
for (1.02 1), but using the Cu/ spin wave velocity. 



3 The Cuu-like modes at (1/2 1/2 0) 

The discussion so far has been of the excitations ob- 
servable near the reciprocal lattice vector (1 1), the 
Cu/ magnetic Bragg peak position. Unlike the spin- wave 
energy, which depends only on the reduced wave vector 
q = Q — G, the neutron scattering intensity from spin 
waves depends also on the reciprocal lattice vector G. 
The neutron scattering intensity is strong near an anti- 
ferromagnetic G, or (w 7r) position, while it is weak near 
a nuclear Bragg position, (0 0). From the spin- wave cal- 
culation, we have found that the Cu//-like modes have 
very large intensity near the Cu// magnetic Bragg po- 
sition, while the Cu/-like modes have vanishingly small 



intensity. Although the Cu//-like modes have large inten- 
sity near (1/2 1/2 L), as illustrated in Fig. |l2|(a), a rather 
complex dispersion relation results due to the presence of 
different magnetic domains (see Fig. 0). Spin- waves from 
domain A are shown as dashed lines, while those from do- 
main B are shown as solid lines in the figure. Therefore, 
one expects to observe three or four peaks within an 1 
meV range around uj = 2.5 meV from neutron scatter- 
ing; this is an extremely difficult task, considering that 
the experimental resolution is about 0.2 ~ 0.3meV in this 
energy range with cold neutrons. 

Representative scans are plotted in Fig. |l2|(c); the data 
have been taken at the SPINS spectrometer at the NCNR 
with collimations of 30'-80'-S-80'-100', and with the fi- 
nal neutron energy fixed at 5 meV. The solid line and the 
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FIG. 12: (a) Spin-wave dispersion of the Cu//-like modes 
along the L-direction near the Cu// antiferromagnetic 
Bragg position (1/2 1/2 L). Since different magnetic 
domains have different reciprocal lattice vector (see Fig. 
|3j), rather complicated dispersion relations result. Spin- 
waves from domain A are shown as dashed lines, while 
those from domain B are shown as solid lines, (b) 
Constant-w scan of the Wi mode near the (1/2 1/2 4) 
position. Solid lines are fits to Eq. (g), and the fitted L 
value of the peak is plotted in part (a) as open diamonds. 
(c) Constant-Q scans of the Cu/j-like modes. 



dashed line for the L = data are the results of fits to two 
spin-wave modes and one spin-wave mode, respectively. 
The two- mode fit is clearly better than the one- mode fit. 
Fitting current data with three spin-wave modes yields 
no meaningful results, since the error bars are larger than 
the separation in peak energies. Therefore, all data have 
been fitted with two spin-wave modes, and the results of 
this fitting are plotted as the solid and open diamonds 
in Fig. [l2](a). The solid circles in Fig. |l^(a) are the fit- 
ting results from coarse resolution measurements using 
higher energy neutrons. A constant-w scan is shown in 
Fig. H|(b) near the (1/2 1/2 4) position at u = 0.8meV. 
The solid line is a fit to Eq. m). Since u>i is an in-plane 
mode, we are able to observe this at large L values. 
From these data, we can establish the following: First, 



the constant-Q scan reveals that there is more than one 
mode in the 2 to 3 meV energy range, roughly coinciding 
with the theoretical prediction. Note that the theoret- 
ical prediction, shown as solid and dashed lines in Fig. 
|l2|(a) , is obtained with parameters determined from pre- 
vious sections, and thus contains no adjustable parame- 
ters. Second, in agreement with the theoretical predic- 
tion, one of these modes is an in-plane mode and the 
other an out-of-plane mode. However, more experiments 
with higher resolution will be valuable in understanding 
the observed spin- wave dispersion. 



D Spin-wave dispersion of Cuj; in the plane 

Because J]j is relatively small, the Cu// zone- 
boundary spin-wave energies are low enough to be ac- 
cessed with thermal neutrons. We have measured spin- 
waves in the afr-plane, along the high-symmetry direc- 
tions. The experiment was conducted at 10K, which is 
well below Tjv,/j, and in the (H K 0) zone; that is, ab- 
plane is in the scattering plane. Both constant-w scans 
and constant-Q scans were carried out. Some typical 
constant-w scans, along the [1 1 0] direction, are shown 
in Fig. |lj|(a). The solid lines are obtained from a least- 
square fit to the cross section, convoluted with the instru- 
mental resolution. Examples of constant-Q scans along 
the zone boundary are shown in Fig. 13(b). 

Our zone boundary data exhibit a double-peak struc- 
ture. The feature at the high energy is the spin-wave, 
while the low-energy feature is a phonon. We have ver- 
ified the different nature of the scattering of the two 
features by measuring the respective peak intensities at 
equivalent reciprocal lattice positions with larger |Q| 



Such measurements at (| 1 0), (| 1 0), (| 1 0), and 
(| 1 0) are shown in Fig. |l3](c): The intensity of the 
low-energy feature increases approximately as ~ |Q| 2 , 
characteristic of phonon scattering, while the intensity of 
the second feature is nearly independent of |Q|. 

Figure |l4J summarizes our results. From the zone 
boundary spin-wave energy of 25meV one can deduce 
Ju rather accurately as Ju = 10.5(5)meV, in excellent 
agreement with the value deduced in Ref. E(J from the 
Cu// susceptibility. The gap energy at the zone cen- 
ter, ~ 3meV, corresponds to the modes found at L = 
in Fig. [12. Away from the 2D zone center, u>\ and 0J2 
from both domains are degenerate and can be approx- 
imated as the excitations of a simple SLQHA with the 
exchange interaction Jjj. The long- wavelength effects 
of the spin-wave interactions can be absorbed into an 
effective anisotropy otf/- With an ~ 0, Eq. (||) can 
be interpreted as resulting from an effective anisotropy 

given by Juotjj = 2Jiai8/{AJjaj + 8) w O.lmeV, or 
eff 



0.01. 



Simple linear SWT with aV — 



'11 



0.01 and J n = 
10.5meV gives the dashed line in Fig. |lj. This is a good 
approximation, except for the dispersion near the zone 



15 



(a)(HHO) 



(b) Zone boundary 



I ' I ' I ' 
"C0=18meV 





FIG. 13: (a) The constant-w scan of Cun spin-waves at 
T = 10K along the high-symmetry direction (H H 0). 
The 2D zone center (w ii) and zone boundary (tt/2 tt/2) 
correspond to Q=(0.5 0.5 0) and (0.75 0.75 0), respec- 
tively, (b) Constant-Q scans along the zone boundary. 
The (1.5 1.0 0) is a local minimum and corresponds to 
(it 0). (c) The \Q\ dependence of each peak in (b) at the 
(1.5 1.0 0) or equivalent positions, which confirms the 
low-energy (20meV) feature as a phonon. 

edge (ir 0). As seen by the continuous line, our data are 
in much better agreement with axecent series expansion 
prediction by Singh and Gelfand.E3 This theory predicts 
a local minimum at the zone boundary position (tt 0), 
lower by about 7% than the value at (| § ). A non- 
zero dispersion along the zone boundary may also result 
from a non-zero next-nearest-neighbor interaction J"™ n , 
within linear SWT. The magnitude of the dispersion be- 
tween (tt 0) and (^ -|) is given by 2SJ^ I n . Considering 
that J a is already of order 10 meV and the next-nearest- 
neighbor distance is large (~ 7.7A), it is unlikely that the 
next nearest neighbor effects contribute strongly to the 
observed zone-boundary energy difference of ~ 2meV in 
Sr 2 Cu 3 4 Cl 2 . 

Therefore, this dispersion can be regarded as a pure 
quantum mechanical effect for the S =f-1/2 nearest neigh- 
bor Heisenberg model. Canali et alxiQ obtained similar 
but smaller zone boundary dispersion in their higher or- 
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FIG. 14: Cu/7 in-plane spin-wave dispersion at T — 10K. 
The dashed line is a fit to the linear SWT with Jn — 
10.5meV and a e // ra 0.01. The solid line is the series 
expansion result (Ref. 52) with the same values for Jn 
and an. 



der SWT. They calculated the correction to Z c up to 
\/S 2 order, and found that the correction is not uni- 
form along the zone boundary, giving ~ 2% dispersion. 
In their spin-rotation-invariant theory, Winterfeldt and 
IhleH also obtained a local minimum at the (7T 0) position 
which is smaller than the energy at the (? ^) position 
by almost 10%. Recenk-quantum Monte Carlo study by 
Syljuasen and R0nnowE3 also gives similar zone bound- 
ary dispersion of 6%, in good agreement with the series 
expansion result and our experimental result. 



V. MAGNETIC CORRELATION LENGTH 

The static structure factor provides valuable informa- 
tion about thermodynamic quantities such as the cor- 
relation length. As discussed in Sec. ||, the necessary 
energy integration can be done automatically, in low- 
dimensional systems, via a 2-axis neutron scattering tech- 
nique. In this section, we present our neutron scattering 
results from such 2-axis measurements, for both the Cu/ 
and Cujj subsystems at temperatures higher than their 
respective Neel temperatures. 



A Cu; system 

The magnetism above Tjv,j of the Cuj system is essen- 
tially the same as that of La2Cu04 or Sr 2 Cu02Cl2. The 
difference in the inter-plane coupling is not important at 
temperatures well above Tnj. The only difference is the 
antifcrromagnetic superexchange J/, which is estimated 
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TABLE III: The Cu-O-Cu superexchange interaction 
and lattice constants of cuprates 



(a) CUj 



(b) Cu T 



La 2 Cu0 4 Sr 2 Cu02Cl 2 Sr2Cu 3 4 Cl2 



a (A) 
Ji (meV) 



3.81 

132(4)° 



3.967 
125(6) * 



3.86 
130(5) 



to be 132(4)mciLfor La2Cu04Jjrom the neutron scatter- 
ing experiment !E\ Greven et alxl extracted 125(6)meV for 
Sr 2 Cu02Cl2 from the two magnon Raman scattering ex- 
periment by Tokura and coworkers^ The Cu-O-Cu su- 
perexchange energies as well as in-plane lattice constants 
in these materials are compared in Table III, where the 



value for Ji in Sr2Cu304Cl2 is extracted from our corre- 
lation length data. 

The cross section for an energy integrating scan across 
the Cuj 2D fluctuations is given by 

I{q2 D )~ / <L>S(Q,u>) (16) 



« [(sin 2 ( f ) )S T (q 2D ) + (sin 2 + 2cos 2 ( f>)S L {q 2D )i 

where q 2 D = — \H — 1|, and S T and S L are the trans- 
verse and longitudinal components of the static fluctu- 
ation. Here <j> is the angle between Q and [0 1]. At 
temperatures well above T^.i only the Heisenberg term 
is relevant in the Hamiltonian. In this regime the system 
is effectively isotropic; that is, S T w S L . 

The neutron scattering data shown in Fig. nB were 
obtained with the incoming neutron energy fixed at 
Ei = 36.4 meV and a collimation sequence of 10'— 13'- 
S-10'. Higher order neutrons were filtered by both PG 
and sapphire filters. The data were fitted to a simple 2D 
Lorentzian convoluted with the instrumental resolution: 

S{Q2d) = -r^-, — 2> ( 17 ) 

q 2 D + k 2 

where the width of the Lorentzian k = £ _1 is equal to 
the inverse correlation length. 

The fitting results for the inverse correlation length, k, 
are shown in Fig. ^q(a) as open diamonds. We have also 
used Ei = 13.7meV neutrons to improve the resolution 
at lower temperatures; these results are shown as open 
circles in the same plot. We alsQ.plnt.the quantum Monte 
Carlo data from several studies .Boo The solid line is the 
renormalized classical (RC) expression of the QNLctM:l£| 
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(18) 

where p s is the spin stiffness constant and v is the spin- 
wave velocity. For the S = 1/2 SLQHA, a recent Monte 
Carlo study by Beard et alE obtains p s /J = 0.1800(5) 
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FIG. 15: Representative energy-integrating two-axis 
scans. The solid lines are fit to simple Lorentzians, 
Eq. (p7|). (a) Scans across the Cu/ 2D magnetic rod 
at (H 0.327); the dashed line shows the instrumental 
resolution, (b) Scans across the Cujj 2D magnetic rod at 
(H H 0.253); the dashed lines indicate the temperature 
independent background. 



and v/Ja — 1.657(2), and we have used these values 
substituted into Eq. ( jig ) to obtain the solid line in 
Fig. p9(a). Although our data have large error bars, 
the general agreement between our experimental results 
and both theoretical results is quite good. We extract 
the value of J/ by comparing our data with quantum 
Monte Carlo results: Jj = 130(5) meV. The experimen- 
tal data deviate from the 2D Heisenberg prediction as 
the temperature approaches Tn,i from above, since the 
system crosses over to 3D Heisenberg behavior due to 
the inter-plane coupling J/.3D- We also show the fitting 
results for the Lorentzian amplitude Sqk 2 — So/t; 2 in 
Fig. Qb). For the QNLcrM,! this quantity is predicted 
to behave as ~ T 2 at low-temperatures, while various 
neutron scattering studiesQa reveal the empirical behav- 
ior So/£ 2 ~constant. Our data are compared with these 
two scaling behaviors in Fig. |l6|(b). The solid line is 
~ T 2 , while the dashed line is a constant; these lines are 
rescaled to fit the data. One should, however, note that 
the low temperature data (below 400K) probably do not 
show true 2D Heisenberg behavior but a crossover to 3D 
Heisenberg behavior, and they should be excluded in the 
comparison. Within experimental error bars, both lines 
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FIG. 16: (a) Inverse magnetic correlation length and (b) 
Lorentzian amplitude (So/£ 2 ) as obtained from fits of the 
Cuj static structure factor to a single 2D Lorentzian. Dif- 
ferent open symbols denote different experimental config- 
urations. In (b), they are scaled to match in the overlap- 
ping temperatuse-^jange. Solid squares in (a) are Monte 
Carlo results pEaEil and the solid lines arc theoretical pre- 
dictions from the QNLcxM with J — 130meV. 



FIG. 17: (a) Inverse magnetic correlation length and (b) 
static structure factor peak amplitude So as obtained 
from fits to a single 2D Lorentzian. Different open sym- 
bols denote different experimental configurations. In (b), 
they are scaled to match in the overlapping tempnmture 
range. Solid squares are Monte Carlo results,tltBo and 



the solid line at T < 40K in (a) is k 
described in the text. 



0.026A- 1 as 



describe our data equally well. 



B Cuh system 

The Cu/7 two-axis cross section is given by 
I{q2 D )= / dwS(Q,w) (19) 



[(2 + sm 2 ( ( f>))S T (q 2D ) + (1 + co S 2 (0))S L (q 2D )}, 



where q 2 D 



{H-\Y + {K 



±) 2 , and S T and S L 



are the transverse and longitudinal components of the 
static fluctuations. S L diverges at the Ising ordering tem- 
perature Tnji- The particular geometric factors result 
from the fact that the Cu/j easy-axis lies within the cop- 
per oxide layers and that there exist two types of domains 
that are equally probable. This gives an almost 3 to 1 
ratio of transverse to longitudinal components when L is 



small, making it difficult to observe longitudinal fluctu- 
ations. In K2N1F4, where the easy-axis is perpendicular 
to the NiF 2 plane, this ratio is close to 1:1, enabling one 
to observe readily the longitudinal (Ising) component of 
the static structure factor.E3 

Our experiment was carried out with Ei = 14.7meV 
and with collimations 20'-40'-S-40'. Representative 
scans are shown in Fig. |u](b). At higher temperatures, 
we used Ei = 30.5mcV in order to ensure that the en- 
ergy integration is done properly, since the characteristic 
energy scale becomes large at these temperatures. We 
could not distinguish the longitudinal component from 
the transverse component; therefore, the solid lines in 
Fig. [l5|(b) are results of fits to a single 2D Lorentzian, Eq. 
(O), convoluted with the experimental resolution. The 
so-obtained correlation lengths versus temperature are 
plotted in Fig. 0(a). Also shown are MontepGpJo results 
for the S = 1/2 nearest-neighbor SLQHA.BESeJ At tem- 
peratures well above Tnji the spin system is effectively 
isotropic, and the correlation length agrees very well with 
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the numerical result. However, even at lower tempera- 
tures the agreement is quite good, since the transverse 
term in the cross section is 3 times larger than the lon- 
gitudinal term. In addition, Ising criticality has a very 
small critical temperature range. As the temperature is 
lowered, the crossover from the 2D Heisenberg to the 2D 
Ising symmetry presumably occurs very close to the tran- 
sition temperature, TV.//, and hence isotropic behavior 
is observed for T > Tmji- We also show the static struc- 
ture factor peak amplitude S in Fig. |lj](b) along with 
the Monte Carlo results from Ref. |69. Similar to the in- 
verse correlation length data, the agreement is quite good 
for all T > T NJI . 

At temperatures below TV,//, the inverse correlation 
length shows a saturation around 0.025A -1 . In thek 
study of the 2D antiferromagnets, Birgeneau et aZ.EJ 
showed that the transverse susceptibility dominates be- 
low the Neel temperature and it can be described via 
spin wave theory. In the presence of an Ising anisotropy, 
a, a typical spin wave dispersion is given as cu q oc 



8 J 2 Z 2 a + v 2 q 2 /2 and the spin wave intensity is pro- 
portional to 1/wg, where v is the spin wave velocity, 
2S^/2Z c Ja. For lj/T <C 1, the population factor is re- 
duced to Uq 1 . Thus the neutron scattering intensity of 
the wave vector dependent susceptibility is 

' ' ' L r , (20) 



Z 2 J 2 a 2 q 2 



1 



which is a Lorentzian with a finite width, k± = ^4 
By substituting a^j ~ 0.01 for the a, we obtain k± rs 
0.026 A -1 . This value, indicated as a solid line in the 
figure, agrees remarkably well with the experimental re- 
sults. 

In Fig. ftty, the correlation length data for T > 43K are 
plotted as a function of inverse temperature on a semi- 
log scale. Also plotted are the correlation length data 
of Sr2Cu02Cl2 taken from Ref. 0. These neutron scat- 
tering data are compared with various theoretical pre- 
dictions. The RC expression for the correlation length 
is plotted as a dot-dashed line. Quantum Monte Cada 
results and high temperature series expansion resultsEa 
are shown as solid and dotted lines, respectively. In a 
recent theoretical study, Cuccoli et alSn treated quan- 
tum fluctuations in a self-consistent Gaussian approxi- 
mation, separately from the classical contribution. This 
purely-quantum self-consistent harmonic approximation 
(PQSCHA) result is plotted as a dashed line in Fig. UH. 
The combined experimental data span almost two orders 
of magnitude in correlation length and show quantitative 
agreement with the Monte Carlo results without any ad- 
justable parameters. At high temperatures, £/a < 10, 
both the series expansion and the PQSCHA, which cor- 
responds to classical scaling, agree with the experimental 
data within error bars. The surprisingly good agreement 
between the neutron scattering data and the renormal- 
ized classical prediction even up to a very high temper- 
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FIG. 18: The logarithm of the reduced magnetic corre- 
lation length £/a versus J/T. The open and filled circles 
are the data for Sr 2 Cu0 2 Cl 2 (Ref. 7) and Sr 2 Cu 3 4 Cl2 
(Gujj), respectively. The RC prediction of the QNLoM is 
plotted in dotrdashed line. Interpolated quantum Monte 
Carlo teaultsQ'Oa and high temperature series expansion 
resultseS are shown in solid andxlptted line, respectively. 
Cuccoli et aUs PQSCHA resultH is plotted as a dashed 
line. 



ature turns out to be a fortuitous one. Beard et al.Q 
pointed out that the renormalized classical scaling sets 
in only at large correlation lengths so that the tempera- 
ture range probed by the neutron scattering experiment 
(T > 0.2 J) is not low enough to see this asymptotic scal- 
ing behavior. However, the deviation is smaller than the 
experimental errors, making it difficult to discern any 
discrepancies from the neutron scattering experiment. 

There are two other recent neutron scattering stud- 
ies on the magnetic correlation length of the S — 
1/2 SLQHA. Birgeneau et alB extended previous work 
on La 2 Cu04 to higher temperature and showed that 
the data are well-described by the Monte Carlo, the 
PQSCHA, and series expansion results within the ex- 
perimental uncertainties. They also showed that there 
is no evidence for a crossover from renormalized classi- 
cal to quantum critical behavior, at-ieast from the cor- 
relation length data. R0nnow et aliB also have carried 
out a study of the correlation length in the monoclinic 
planar antiferromagnet copper formate tetra-deuterate 
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FIG. 19: The Lorentzian amplitude of the structure fac- 
tor, So(t;/a)~ 2 , is plotted as a function of T/J. The open 
and filled circles are the data for La2Cu04 (Ref. 8) and 
Sr 2 Cu 3 04Cl2 (Cuj/), respectively. The RC prediction 
of the QNLcrM is plotted as a solid line. Our quantum 
Monte Carlo results are also shown in open squares. 



(CFTD). They obtained essentially similar results to 
those shown here for the Cu/j system, agreeing with 
the Monte Carlo data up to a very high temperature 
(T w 1.25 J). They were able to extend their measure- 
ment to such a high temperature by employing a special 
technique involving filtering out the elastic part of the 
signal, thus reducing the incoherent background. How- 
ever, the analysis depends sensitively on the theoretical 
model, especially on the scaling of the characteristic en- 
ergy scale r q= o, which still needs further investigation. 
Besides, its low-symmetry crystal structure and relatively 
large Dzialoshinsky-Moriya interaction make CFTD a 
less ideal S = 1/2 SLQHA than Sr 2 Cu304Cl 2 (Cu//). In 
fact, the combined Sr 2 Cu0 2 Cl 2 -Sr 2 Cu304Cl 2 (Cu/_r) sys- 
tem forms an ideal model 5 = 1/2 SLQHA over a large 
temperature range 0.2 < T/J < 0.75. 

In Fig. H9, the Lorentzian amplitude of the struc- 
ture factor, So(£/a)~ 2 , is plotted as a function of T / 
Our data and the La2Cu04 data of Birgeneau et al 
are scaled to match the Monte Carlo results, which 
are plotted in absolute units without any free parame- 
ter. The RC prediction is also plotted as a solid line: 

|£ = A2-kM 2 s (2^), with A s=1/2 = 3.2 from the se- 
ries expansion studyO The first thing to note in our data 
is the disappearance of the divergence at Tjyjj, which 
implies that the divergence in So is absorbed by the £ 2 
term, or equivalently r\ = 0, as predicted for the 2D 
Heisenberg model. On the other hand, the critical ex- 
ponent r\ for the 2D Ising model is exactly known to be 
7] = 1/4; thus So£~ 2 should show a weak divergence of 
£_i/4 _ ^ _ T )-i/4 ; which is not b S erved. This is 



not surprising, since the finite Q-resolution prevents us 
from observing even the strong divergence of £ in the first 
place. 

What is surprising, however, is the discrepancy ob- 
served at high temperatures between the two experimen- 
tal sets of data. Unlike La 2 Cu04, which shows constant 
So£ over the observed temperature range, the Cu// 
system shows some temperature dependence. Specifi- 
cally, So£ -2 follows the Monte Carlo data closely for 
T > 0.4 J and deviates significantly from the RC predic- 
tion at high temperatures. S £~ 2 for Sr 2 Cu0 2 Cl 2 shows 
behavior similar to that of La2Cu04, while CFTD shows 
behavior similar to that of Sr2Cu304Cl2. Considering 
that the CETD also has rather small exchange coupling of 
~ 6.3meVjI3 one can speculate that this discrepancy may 
be due to high energy fluctuations at high temperatures, 
which are not integrated in the experiments on systems 
with a large exchange interaction, such as La2Cu04. An 
experimental study on dynamic critical behavior of the 
2DQHA is necessary to address this problem further. 

Recently, in their study of the S = 5/2 SLQHA, 
Rb 2 MnF4, Leheny et alE3 showed that the amplitude 
So£ -2 crosses over from constant behavior at high tem- 
perature to the ~ T 2 behavior at low temperature. The 
crossover temperature was estimated to be around ~ 
7.4 J, which would correspond to ~ 0.6 J for the S = 1/2 
case, after the temperature scaling of S(S + 1) is taken 
into account. One might be able to regard the high tem- 
perature data for T/J > 0.6 of So£~ 2 as constant, thus 
showing similar crossover behavior, although it is difficult 
to make any strong statement with the current data. 



VI. DISCUSSION 

From our study, we have been able to determine two 
very important superexchange interactions: one is the 
"edge-sharing" Cuj-O-Cuj/ exchange interaction, and 
the other is the Cujj-Cujj interaction, which corre- 
sponds to the second nearest neighbor interaction in the 
Cu/ square lattice. First, we estimate the isotropic Cuj- 
Cu// interaction as Ji-u ~ — lOmeV. This edge shar- 
ing superexchange interaction is crucial in understanding 
spin ladder materials, SrCu 2 03 and Sri4Cu 2 404i, as well 
as other ID spin systems, such as SrCu0 2 (S = 1/2 zig- 
zag chain) and CaCu 2 03 (buckled ladder). It has been 
assumed that the edge sharing interactions, which hap- 
pen to be frustrated in all these materials, are small, and 
that they therefore can be ignored in data analysis. How- 
ever, as we have seen in Sr 2 Cu304Cl 2 , quantum effects 
may be important in understanding the low temperature 
properties of these systems. 

The second nearest neighbor interaction in the copper 
oxide plane of the high temperature superconductors has 
the same superexchange path as the Cu//-Cu// inter- 
action; namely, the Cu-O-O-Cu path. Of course, the 
Cujj-Cujj interaction has additional contributions from 
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the path Cu-O-Cu-O-Cu; however, one would expect 
this contribution to be small in magnitude; it should also 
be ferromagnetic. Thus, we can assume that the second 
nearest neighbor coupling in the copper oxide plane must 
be close to the Ju value: ^lOmcV. This value is used to 
fit the ARPES data m,thc framework of the t-t'-t"-J 
model by Kim et aZ.Jl3 which shows good agreement. 

One expects interesting physics to arise from dop- 
ing this system with either charge carriers or non- 
magnetic impurities. One difficulty of studying doped 
Sr 2 Cu 3 04Cl2 is that it is extremely difficult to dope this 
system with any impurities. Many attempts to dope this 
system with impurities such as Zn, Mg, K, Y, etc. have 
failed. In fact, there are only two successfully doped 
cogper oxy-hallides: Sr2Cu02F2+<5 by AUMamouri et 
alB and (Ca,Na) 2 Cu0 2 Cl 2 by Hiroi et al.B Both com- 
pounds are superconducting and are synthesized at high- 
pressure. 

In his study of a frustrated vector antifcrromagnet, 
where the next nearest neighbor coupling is much greater 
than the nearest neighbor coupling, HenleyO has shown 
that the disorder introduced by dilution favors anti- 
collinear ordering. Since quantum fluctuations prefer a 
collinearly ordered ground state, these two types of dis- 
order compete with each other and produce a rich phase 
diagram as a function of temperature and dilution. How- 
ever, diluted Sr2Cu304Cl2 is a little different in that the 
relevant coupling ratio (Ji-n/ Ju) is not small, so that 
the simple perturbation expansion used by Henley is no 
longer applicable. Nevertheless we expect a dramatic 
change in the ground state of diluted Sr 2 Cu304Ci2; for 
example, a helical order, a spin glass, or even a disordered 
ground state might occur as a result of dilution. 



VII. CONCLUSIONS 



tion between these subsystems is frustrated. Some modes 
(i. e. u>s) would have zero energy in the absence of quan- 
tum fluctuations. Other modes (i. e. 102 and L04) show 
remarkable effects of quantum fluctuations. In all cases, 
these dramatic shifts are in quantitative agreement with 
theoretical calculations. 

(4) Our measurement of the spin-wave dispersion al- 
lows precise determination of several exchange inter- 
actions, including the interplanar Cuj-Cur interaction 
(Ji,3D = 0.14(2) mcV), the in-plane Cu//-Cu// inter- 
action (Ju = 10.5(5) meV), and the in-plane Cuj-Cu// 
interaction (\Ji-u\ = 10(2) meV). 

(5) We have made precise tests of spin-wave interac- 
tions at the zone boundary which support recent theo- 
retical calculations. This test is particularly convincing 
for this spin 1/2 system, where these effects are too large 
to be attributed to further-than-nearest neighbor inter- 
actions. 

(6) The instantaneous spin-spin correlation length, £, 
of 5 = 1/2 SLQHA over a wide temperature range has 
also been obtained from our neutron scattering experi- 
ments. Our measured values of £ are in good agreement 
with recent calculations based on the quantum nonlinear 
a model and on quantum Monte Carlo simulations. 
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We have presented results from our neutron scattering 
experiments on S^CuaO-jC^, and discussed its magnetic 
properties as well as the novel quantum phenomena as- 
sociated with order from disorder. In what follows, we 
briefly summarize our main results. 

(1) Our elastic neutron diffraction data confirm the. 
magnetic structure obtained from a previous analysis!^ 
of static properties in a magnetic field. 

(2) We show that Sr2Cu304Ci2 is a unique system hav- 
ing two independent phase transitions. By analyzing the 
intensities of the magnetic Bragg reflections we obtain 
the critical exponents for the order parameter for the 
Cu/ transition at T NJ = 386(2)K fa = 0.28(3) and for 
the Cu// transition at T NJI = 39.6(4)K /3 n = 0.13(1). 
The Cu/ transition is thought to be that of a 3D XY 
model, whereas the Cu// transition is identified as a 2D 
Ising transition. 

(3) The dramatic variation in the mode energies as the 
Cu// subsystem orders is very clear evidence of quantum 
fluctuations because on the mean field level the intcrac- 
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